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I. INTRODUCTION 

The magnetic properties of materials are usually de- 
scribed in the frame of system models, such as Heisenberg 
or Ising models where rigorous results, or suitable mean 
field approximations are available in the thermodynami- 
cal limit. On the other side, modern applications require 
to deal with nano-sizcd magnetic materials, whose intrin- 
sic features lead, from one side to the emergence of quan- 
tum phenomena^] , and to the other to the question of 
applicability of statistical mechanics. Indeed, few parti- 
cle systems do not usually fit in the class of systems where 
the powerful tools of statistical mechanics can be applied 
at glance. In particular, an exhaustive theory able to fill 
the gap between the description of 2 and 10 23 interact- 
ing particles is still missing. Moreover, also important 
well-established thermodynamical concepts as the tem- 
perature, become questionable at the nano-scaleQ- 

In a similar way, long-range interacting systems be- 
long, since long, in the class where standard statistical 
mechanics cannot be applied tout court. Indeed, they dis- 
play a number of bizarre behaviors, to quote but a few, 
ensemble inequivalence negative specific heat, temper- 
ature jumps and long-time relaxation (quasi-stationary 
states) 0. Therefore, from this point of view, few-body 
short-range interacting systems share some similarities 
with many-body long-range ones. 

Within such a scenario, and thanks to the modern com- 
puter capabilities, it is quite natural take a different point 
of view, starting investigations directly from the dynam- 
ics, classical and quantum as well. It was in this spirit 
that, few years ago, a topological non-connection of the 
phase space was discovered £| in a class of anisotropic 
spin systems. This was initially called, for historical 
reasons @, breaking of ergodicity, meaning with that a 
trivial consequence, namely that the system can not be 
ergodic (the phase space is exactly decomposable in two 
unconnected parts) [3, even if we prefer here to call it 
Topological Non-connectivity Threshold (TNT). This re- 
sult, was found first numerically and later analytically, 
in a class of models, the anisotropic Heisenberg models, 
where important and rigorous results have been obtained 



during the last century, in the thermodynamical limit 
only. 

Quantum effects in such small magnetic systems can 
not be neglected, in principle, even if the usual viewpoint 
is to consider magnetic domains as quantum objects 
with huge spin number. Still, what we have in mind now 
and in our future plans, is to show the relevance of the 
TNT with respect to the complicated transition between 
the classical and quantum world. For instance, it is well 
known that quantum particles can tunnel across potential 
barriers at variance with the classical ones. What is less 
obvious is that a macroscopic variable, such as the mag- 
netization, can do the same. This phenomenon, known 
as Macroscopic Quantum Tunneling, well described in 
n| is an important step in the so-called Leggett program 
; S,J for a better comprehension of the classical-quantum 
transition. Thus, after a brief description of the quantum 
analogue of the classical TNT, we show its relevance in 
single-spin models used in micromagnetism, featuring the 
TNT as a perturbative threshold. 



II. THE QUANTUM TOPOLOGICAL 
NON-CONNECTIVITY THRESHOLD 

The results found in the classical model [jj IqL lioll has 
been considered in the semiclassical regime in [l II . 

Here, we consider a system of N particles of spin 1 = 1, 
described by the following Hamiltonian: 



N 1 JV 



(1) 



=1 j^ii 



=1 &i 



where — 1 < 77 < 1 is the anisotropy constant. Quantiza- 
tion of the Hamiltonian follows the standard rules. (Let 
us remember that, according to the correspondence prin- 
ciple, the classical limit is recovered as I — > 00). As in 
the classical case we fix the modulus of the spins to one. 
This can be achieved with an appropriate rescaling of the 
Planck constant, h -> h/\Si\ = 1/^1(1 + 1). With this 
choice, in the classical limit, I — ► 00 (h — > 0), the spin 
modulus remains equal to 1. We will also limit our anal- 
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FIG. 1: Energy splittings S versus E for the many-spin case. 
Eigenvalues and splittings are compared for the case N — 20, 
rj = 1 and 1 = 1. Circles : mean field approximation; squares 
: full Hamiltonian; crosses : perturbative result. Eigenvalues 
are arranged in regular block, according to different L values 
as indicated in the figure. 
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FIG. 2: Probability p = \(E\Emf)\ 2 that an eigenstate \E), 
with specific energy E/N, occupies an eigenstate \Emf), with 
specific energy Emf/N. Parameters are : a) N = 100, I = 1 
;b) N = 6l = 4. 



ysis in the subspace of all possible completely symmetric 
states (bosonic symmetry) . 

In |ll| it was shown that the magnetization along the 
easy axis, at variance with the classical case, can change 
its sign below the TNT through Macroscopic Quantum 
Tunneling. This leads to the problem of a significant def- 
inition of the quantum TNT. In the semiclassical limit 
(large I) a quantum signature of the classical TNT can 
be found pd| in the spectral properties of the system lead- 
ing to a proper definition of the quantum disconnection 
threshold, E^ nt , with the correct classical limit. Below 
Ef nt the spectrum is characterized by the presence of 
quasi degenerate doublets, whose energy difference, S, 
increases exponentially up to Ef nt , and saturates above 

a tnf 

On the other side here, we focus on the hard quantum 
regime (/ = 1). The energy spectrum, still presents dou- 
blets and an approximate exponential dependence of 5 
with the energy. Nevertheless, it is evident that, at vari- 
ance with the semiclassical case, they change regularly 
by many order of magnitude in small energy bins, see 

pig. ni 

In order to understand the origin of this regularities, 
it is useful to rewrite JTJ as H — Hm f + H\ , where 

T) 1 

Hmf = -rh x - -rh y (2) 

Bi = -E^T-lE^f) 2 - (3) 

and rh XtVtZ = ^ i S^' y ' z . While the first term (mean 
field) is integrable in the classical limit, the latter is re- 
sponsible for the non integrability of the system. Let 
us also consider the eigenvectors of Hmf, \Emf), and 
expand the eigenvectors of H, \E) over them. In other 
words we consider the probability, po = | (Em f | E) | 2 , that 



a given eigenvectors of H occupies a given eigenvector 
of Hmf- As one can see, see Fig. the eigenvectors 
of the full Hamiltonian are almost completely localized 
on the eigenvectors of the mean field Hamiltonian. over 
the whole energy range. Actually in the low energy re- 
gion the eigenvectors occupies just one eigenstate of the 
mean field Hamiltonian with probability greater the 0.9, 
while all the other states are occupied with probability 
smaller the 0.01. The same does not happen in the large 
I case (Fig. |2Jd). Therefore, the non- integrable part is 
negligible with respect to the mean field J2J- The ques- 
tion of the quantum integrability of chaotic Hamiltoni- 
ans for bosons with I = 1 has been recently posed in 
[l2|. Shortly, quantum integrability should be induced, 
for I = 1, by the strong correlations between Hamiltonian 
matrix elements. From the above analysis it follows that 
we can use the mean field Hamiltonian to study the total 
Hamiltonian in the hard quantum regime. 

Here, we present the results of a high order pertur- 
bative calculation of the eigenvalues of Hmf- Since 
[HmFtIti 2 ] = 0, it is sufficient to consider the eigenvalues 
of to 2 . They are given by the possible values of the to- 
tal magnetic moment which can be obtained combining 
./V particles of spin I, and are determined by the quan- 
tum numbers: Nl,Nl ~ 1,....,0. From these values we 
should exclude those which cannot be combined to give 
completely symmetric states, if one is interested in the 
bosonic case (even if the present approach is indepen- 
dent from the statistics) . We can consider each subspace 
with different m 2 separately. In this way the many-spin 
Hamiltonian Hmf, is equivalent to a set of single spin 
systems, described by the same Hamiltonian. Note that 
I is the magnitude of the spin of the many-spin problem. 
Thus in the following we will first consider single spin 
models, and then we will come back to our many-spin 
problem. 



III. SINGLE SPIN MODEL 

Let us consider a single spin of magnitude L, z— com- 
ponent L z , with \L Z \ < L. For simplicity we will rewrite 
the mean field Hamiltonian as follows: 

Hmf = - {-ml + f]m 2 x ) (4) 

Eq. (0} can be reduced to after a rotation of it around 
the x axis which carries y in z and z in — y, and does not 
affect the physics of the problem. 

Single-spin-model have an interest in themselves, be- 
sides the fact that their analysis will allow us to compute 
the energy levels of our many-spin mean field Hamilto- 
nian. In recent years growing interest arose in micromag- 
netic particles |lill3|. such as ferromagnetic domains and 
magnetic macro-molecules such as Mnn and Fes EH- 
The research interest in these systems is mainly due to 
the possibility to reveal quantum effects in the macro- 
scopic domain, such as the Macroscopic Quantum Tun- 
neling (MQT) of the magnetic moment, and the even 
more interesting phenomenon of Macroscopic Quantum 
Coherence (MQC). While in the former case (MQT) the 
total magnetization of a microscopic particle flip even 
if classically this would be forbidden by the presence of 
an effective energy barrier, in the latter case (MQC) the 
magnetization oscillates between opposite magnetization 
states in a coherent way. This phenomenon, if revealed, 
would unambiguously indicate the presence of Quantum 
Interference of Macroscopic Distinct States || . At suffi- 
ciently low energy these systems can be modeled by phe- 
nomenological single-spin Hamiltonians, where the single 
spin describes the total magnetic moment of the system. 
Splittings of the eigenvalues of the single spin Hamiltoni- 
ans are simply related to the frequencies of MQC (or the 
tunneling rates of MQT) 1]. For this reason much effort 
has been devoted in these years to compute such split- 
tings. Usually, semiclassical methods are employed, such 
as WKB and imaginary time path integrals to quote but 
a few ^| . Also perturbation theory can be success- 
fully applied in this kind of problem taking into account 
high order terms |17| . Indeed it is possible to compute 
explicitly the first non-zero perturbative contribution to 
the splittings even if this is an high-order contribution. 

In Appendix [V] we show the basic ideas of the high 
perturbative order approach, and we derive an analyti- 
cal expression for the eigenvalues and the corresponding 
splittings (in 17] no explicit derivation was given). 

Numerical eigenvalues and their splittings for Hamil- 
tonian Q and a single spin of magnitude L, have been 
compared with our results in Fig. |3 where we show the 
splittings S as a function of the energy E. In the same 
figure we can see that while the high perturbative order 
approach gives a very good estimate for E < E^ nt it fails 
completely above the threshold. An upper bound to the 
energy at which our approach fails can be given evalu- 
ating the quantum correction to E tn t for the single spin 
Hamiltonian, see [13 (indicated in Fig. as a vertical 
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FIG. 3: 5 shifts versus the energy E. Eigenvalues and split- 
tings obtained numerically (open circles) and via perturbation 
theory (crosses) for a single spin L = 20. Also shown as verti- 
cal dashed line the rough estimates for the validity of the our 
approach. Here is Ti = 1, 77 = —0.5. 



dashed line). 



IV. MANY SPIN HAMILTONIAN 

We now compare our analytical results with the eigen- 
values of the mean field Hamiltonian Hmf considered as 
a many spin Hamiltonian and with the eigenvalues of the 
complete Hamiltonian in the hard quantum regime 
We can achieve this by considering the possible eigenval- 
ues of m 2 , which are obtained when an ensemble of N 
particles of spin I is considered. Note that from the set 
of possible eigenvalues of m 2 we have to exclude those 
that are not compatible with the symmetrization pos- 
tulate (for the bosonic case). For each possible value of 
m 2 we apply our perturbative approach to the correspon- 
dent single spin problem. Then, putting all together, we 
obtain the results for the many-spin Hamiltonian. 

In Fig. n we plot the splittings versus the specific en- 
ergy for the mean field Hamiltonian (circles), for the 
full Hamiltonian Q (squares), and the perturbative re- 
sults (crosses). As one can see we can give an accurate 
good approximation in the low energy region of the spec- 
trum. Deviations obviously appear when the perturba- 
tive approach is compared with the splittings of the full 
Hamiltonian, even if, in this case, the perturbative ap- 
proach still give a correct order of magnitude estimate 
of eigenvalues and splitting. From Fig. we can also see 
how the regular features of S in the hard quantum regime 
are related to the quantum numbers of the total angular 
momentum. 

One could ask if the same perturbative approach works 
in the semiclassical regime for the total Hamiltonian J5J. 
Perturbation theory cannot work in the whole energy re- 
gion, even if we may expect to give an approximate de- 
scription for low-energy eigenvalues. For instance, ap- 
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plying perturbation theory, e.g. Eq. (|14fl . for the energy 
separation between the ground state and the first excited 
state, we get: Sg s = 0.25h 2 /6~ Nl . This expression works 
works fairly well[l3, even for large I and it reproduces 
the main features of the dependence of the ground state 
splitting, namely the exponential dependence on I, and 
N as well. 

We have shown that, in the hard quantum regime, it is 
possible to compute perturbatively the splittings of the 
doublets characterizing the spectrum below the quantum 
TNT. This is due to the nearness of the full Hamiltonian 
with the mean field Hamiltonian. The quantum TNT 
can be also considered a perturbative threshold, since it 
gives an estimate of the energy at which the perturbative 
approach fails. Finally, we point out that this threshold 
indicates an energy range which is not negligible with re- 
spect to the total energy range for long range interacting 
systems in which Macroscopic Quantum Phenomena can 
be studied. 



V. APPENDIX 

In this section we present the results of a high pertur- 
bative order calculation of the eigenvalues of the single 
spin Hamiltonian We will show that in order to 

split the double degenerate levels of the — m 2 term of (@J 



characterized by the quantum number Iq 



the first 



non-zero contribution is at the ^o-th perturbative order. 
This also give a qualitative explanation of the well known 
exponential dependence of the splitting magnitudes on 
the energy. 

Let us consider a single spin of magnitude L, z- 
component L z with \L Z \ < L and \L,L Z ) as basis states. 
Hamiltonian |0J can be written as H = Hq + V, where : 



1 



Hq = —m 2 + —r](rh + m + m m + ) 



V — —f](m + rh + + rh rh ) 



(5) 



and m ± = rh x ± irh y . Since Hq is diagonal in the basis 

\L,L Z ), the unperturbed energy Eq = (L,lo\H\L,lo) are 
given by: 

E (L, lo) = -nHl + \h 2 V [L(L + 1) - I 2 )] (6) 

Each unperturbed energy level turns out to be doubly de- 
generate, with eigensubspaces spanned by \L,±l ). The 
first non-zero contribution to the splitting of a degener- 
ate pair occurs at the lo-th order of perturbation theory. 
In order to show that let us define the n— th order per- 
turbation operator [l9| : 



ra-l 



Eq — H 



-V 



(7) 



where = 1- J2e> ^e \ e 'o)(E'o\ is the projector out 
of the considered degenerate subspace. The right linear 



combination of the unperturbed basis vectors \L, ±Zo) (to 
which eigenstates of H tend when V is negligible) can be 
found by diagonalizing the following matrix: 



V 



>(«) 
++ 

00 



V 



,00 

00 



(8) 



where V { ™) = (L, slolV^^L, s'l Q ), s, s' = ±1 and n is the 
minimum order giving rise to two different eigenvalues of 
the matrix ©. A tt~ rotation around the a;— axis leaves V 
unchanged since \L, lo) —> \L, —Iq), and rh — ► rh T . Then 



,(n) 
++ 



,(») 



•pO 1 ) 



unperturbed basis vectors is: 



T 5 !"^ and the right combination of 



±m ) = (\L,l )±\L,-l ))/V2. 



(9) 



Eigenvalues undergo a shift given by : ±V { +1. The 
generic n-th order energy shift, induced by the per- 
turbation is given by: = (±m |P (n) | ±m ). While 
the degeneracy can be removed only by non-zero off- 
diagonal elements, an overall energy shift D can be in- 
duced by non-zero diagonal elements too. 

In order to compute D and S the action of fi^ on 
the basis states \L, ±Zo) should be evaluated. If n = 1 
then = V . In this case the diagonal elements of the 
matrix (JSJ are zero since V can only change — > ± 
2. Off-diagonal elements (L, —Iq\V\L, +Iq) are different 
from zero only when V brings \L, Iq) into \L, — Iq). This 
can happen only if lo = 1. 

En = 2 then V {2) = V((j>/(E - H ))V. Since 



(Eq — Hq) 



\L,l) 



1 - SUa 



E (L,lo) - E (L,l) 



\L,l), (10) 



in order to understand the action of we have to 
apply V twice. Let's consider the diagonal elements: 
Can we take \L,l ) m itself \L,l ), using V twice? Yes: 
VV\L,l ) -> V(\L,l -2) + \L,l + 2)) -> \L,l ) + \L,l - 
4) + \L, lo + 4) + \L, Iq), where the coefficients in front of 
the states have been omitted. Bracketing the final states 
thus obtained with \L, Iq), only the first and the last re- 
main. Then there are two "ways" in which the operator 
can take \L, lo) in itself: if Iq > 1 by the follow- 
ing chain rule : \L, Iq) — * \L,lo — 2) — > \L,Iq), while if 
lo < L - 1 by \L, l ) -> \L, l + 2)^ \L, l ). 

It is now easy to compute the first non zero contribu- 
tions to the overall shift. From 10 we have: 



D=(L,±l \V 



Eo — Ho 



■V)\L,±lo), (11) 



The only contributions D± different from zero, coming 
from the two ways described above, are: 



o (n^ f^M£(L,k-i) t (12) 



5 



respectively for Iq > 1 (£>_) and < L + 1 (£>+). Here, 
we defined / ± (L^ ) = L(L+ 1) - / (i ± 1) and A£± = 
E {L,l )-E {L,l ±2). Thus, L> = £>_+£>+ is the first 
non-zero overall energy shift. 

Let us now consider the off-diagonal matrix elements. 
It is possible to go from \L, Iq) to \L, —lo) using V twice 
only when Iq ~ 2. In this case there is one only way: 
\L,l ) -> |i,Z - 2) -> \L,l -4) = |i, -Z ), the last 
being true only for Iq — 2. It is then clear why the first 
non-zero operator which splits the doublet characterized 
by L, l is the l -th order. From © we have: 



5 = 2(L,l \V 



Eq — Hq 



V 



\L,-l ) (13) 



From Eq. (|13|l there is only one way to connect \L, Iq) 

with \L, —Iq), namely \L, l ) — > \L, Iq — 2) — > \L, —Iq + 

2) — > \L, —Iq). After some algebra, one has: 



5 = 



(14) 



To summarize, for any given degenerate doublet we can 
calculate the overall shift D, Eas.(jl2[l and the splitting 
6, Eq.<HU>. 
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